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Abstract 

We investigate the spin-dependent dynamical response of a quantum ring with a spin-orbit 
interaction upon the appUcation of Unearly polarized, picosecond, asymmetric electromagnetic 
pulses. The oscillations of the generated dipole moment are sensitive to the parity of the occupation 
number in the ring and to the strength of the spin-orbit coupling. It is shown how the associated 
emission spectrum can be controlled via the pulse strength or a gate voltage. In addition, we inspect 
how a static magnetic flux can modify the non-equilibrium dynamics. In presence of the spin-orbit 
interaction and for a paramagnetic ring, the applied pulse results in a spin-split, non-equilibrium 
local charge density. The resulting temporal spin polarization is directed perpendicular to the 
light-pulse polarization axis and oscillates periodically with the frequency of the spin-split charge 
density. The spin-averaged non-equilibrium charge density possesses a left-right symmetry with 
respect to the pulse polarization axis. The calculations presented here are applicable to nano-meter 
rings fabricated in heteroj notions of III-V and II- VI semiconductors containing several hundreds 
electrons. 

PACS numbers: 78.67.-n, 71.70.Ej, 42.65.Re, 72.25.Fe 
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I. INTRODUCTION 



Advances in nanotechnology opened the way for the synthesis of artificial nanostructures 
with sizes smaller than the phase coherence length of the charge carriersi. The electronic 
properties of these systems are dominated by quantum effects and interferences-. Particulary 
interesting are ring structures which served as a paradigm for the demonstration of various 
aspects of quantum mechanics^. Currently available phase-coherent rings vary in a wide 
range in size and particle density^'^'^'^i^'^. On the theoretical side, various features of the 
equilibrium properties of quantum rings are well understood and documented^. Recently 
no-equilibrium dynamics triggered by external time-dependent electromagnetic fields has 
been the subject of xQsea.ich^^^^^^^^^^'^^'^^'^^^. In particular it has been shown that 
irradiations with picosecond, time-asymmetric, low-intensity light fields generates charge 
polarization and charge currents in a qualitatively different manner than in the case of 
applied harmonic laser fields. Currently, asymmetric pulses are producible with a duration 
from few hundreds femtoseconds up to nanosecond82ii22ii22ii2ii2^. The optical cycle of the 
electric field of the asymmetric pulse consists of a short half cycle followed by a much longer 
and weaker half cycle of an opposite polarity. Hence, under certain conditions, the external 
field acts as a unipolar pulse and therefore it is referred to as an half-cycle pulse (HCP). 

In this study we focus on quantum rings as those fabricated out of a dimensional electron 
gas formed between heterojuctions of III-V and II- VI semiconductors. Spin-orbit interaction 
(SOI) is crucial in these materials. The influence of the SOI on the equilibrium properties of 
these rings have already been studied^^i^i (for more recent works we refer to22ii2Ei2Si2ii2^i2^) . 
In this work, we shall consider the spin- dependent non- equilibrium dynamic of the ring 
with SOI driven by HCP's and in the presence of a magnetic flux, a problem which, to 
our knowledge, has not been addressed so far. As shown below the applied pulse triggers 
an oscillating charge polarization with frequencies dependent on the number of particles, 
the strength of SOI, the intensity of the light and the applied static magnetic flux. The 
energy scale is set by field-free eigenfrequency of the ring. Furthermore, it is shown that 
even though the light does not couple directly to the spin (at the intensities considered here) 
the presence of SOI leads to a temporal spin-splitting in the non-equilibrium local charge 
density. The resulting non-equilibrium, local spin polarization is perpendicular to the light- 
pulse polarization axis. It oscillates with the same frequency of the spin-split charge density 
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and thus its time average vanishes. Experimentally, the induced polarization is measurable 
by detecting the associated radiation emission. In fact, the driven rings can be utilized as 
a source for harmonic generation. As shown below the power spectrum is, to some extent 
tunable by an external static field that controls the strength of the Rashba SO coupling. 

This work is organized as follows, at first we shall derive the Hamiltonian of the ring 
with SOI coupled to the HCP's field and in the presence of a static magnetic flux. In Sec. II. 
we discuss the initial carriers' wave functions, energies with and without SOI. In Sec. III. 
we consider the time-dependent dynamics for the system when applying the pulse field. In 
section IV, we present our calculations for the dipole moment of the ring. Detailed numerical 
calculations and discussions are contained in Sec. V. 

II. QUANTUM RINGS WITH SPIN ORBIT INTERACTION 
A. Hamiltonian 

We shall consider a quantum ring (QR) with a spin orbit interaction subjected to a 
time-dependent electromagnetic field. In a minimal coupling scheme the single-particle 
Hamiltonian^^ reads 

Ho = ^ + Vir) - e$ + ^(a X n), + /i^B ■ a, (1) 

where 11 = p + eA, p is momentum operator, e is the charge of the carrier, and A is the 
vector potential of the external electromagnetic (EM) field. The term V^(r) in Eq. ([1]) is the 
potential confining the particles to the QR; the third term in Eq. ([1]) is the scalar potential 
$ of the EM field, the fourth term is the Rashba SOI with the coupling constant an; the 
components of a are Pauli matrices. The last term in Eq. ([T]) is the Zeeman term describing 
the coupling between the electrons' magnetic moment hb and the magnetic component of 
the EM field. 

The electric and magnetic fields E(r, t) = — V$(r,t) — dA{r,t) /dt, and B(r,t) = V x 
A(r, t) are invariant under the local gauge transformations'^^''^^ ^'i^, t) = $(r, t) — and 
A'(r,t) = A(r, t) + Vx(r, Introducing a unitary operator R = exp(— zex(r, i)/^) we find 

= ^[axn'U n' = p + eA'. 
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The transformed Hamiltonian reads 

^' = ^ + ^(r) - e$' + ^(a x U% + ^(H' x A) ■ a. (2) 

Within the Coulomb (or radiation) gauge, i.e. $ = and V • A = we obtain 

H' - [P + ^^f^!^^^^^ + Vir) + e| + ^ X (p + e(A + Vx))]. 

■'^^[(p + e(A + Vx))xA]-a. (3) 



h 

In what follows we employ a plane-wave vector potential A = Aqc*'^'^ '""'^*^ + c.c. and a gauge 
function x(r, t) = — A(t) ■ r, where k and u are the wave vector and the frequency of the 
EM field. Furthermore we note that in our case the light propagates perpendicular to the 
plane of the ring. The thickness d of the ring will be on the order of nanometers. Thus, in 
the present study 1 ^ k ■ r, r ^ c? and the dipole approximation is justified, even though 
the radius a of the ring could be in the micrometer range. With Vx(r, t) = — A(t), and 
= _r . ^ = -r ■ E(t) we find 

H' = Hsoi + Mt), (4) 

where 

^soi = ^ + ^(r) + ^(a X p)„ (5) 
is the Hamiltonian of a quantum ring with spin orbit interaction, and 

Hi{t) = -er-E{t)+fXBB{t)-a. (6) 

Switching over to cylindrical coordinates and integrating out the radial dependence, Hsoi 
attains the fornt^^i^i^i^i^ 

^soi = o * 2 y'^^f + IT! [(TxCOSLp + aySm.Lp)[id^ + —) 

where (9^ = 0o = ^/e is the unit of flux, (j) = Bna'^ is the magnetic flux threading the 
ring, a is the radius of the ring, and ub = 2jiBB /h. In Eq. ([7]) we added a static magnetic 
field in addition to the applied light field. 
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III. SPIN-DEPENDENT CHARGE POLARIZATION INDUCED BY A SINGLE 
EM PULSE 



A. Ground-state wavefunctions and the spectrum of the ring 



The energy spectrum of a QR with SOI has been discussed in several works2^i^i2ii2^. 



Consistent with Ref. 



28| we find for the angular single particle wave functions "^^i^) (^^^ 



index T refers to the transpose) 

where S and n denote the spin and the integer angular quantum numbers and 

are spinors in the angle-dependent local frame and 

= cos(7/2),6T = sin(7/2), 
= — sin(7/2), fo-*- = cos(7/2). 

In absence of the static magnetic field, the angle 7 is given by (cf. Fig. 1) 

tan 7 = —Qr = —ujr/uq, hjjji = 2aR/a, hwQ = h'^/{m*a^). 

The local spin orientations is inferred from the relations 

h 



(8) 



(9) 



(10) 



5(r)i 



(sin 7 cos ifCx + sin 7 sin cpCy + cos 762) , 



and 



S{r)i = — [sin(7r — 7) cos(7r + ip)ex + sin(7r — 7) sin(7r + ip)ey + cos(7r — 7)62 



(11) 



;i2) 



Thus 7 is the angle illustrated in Fig. 1(a). The limit Qr 00, i.e. for a very strong SOI 
coupling, 7 —7t/2 which corresponds to the plane of the ring. The eigenenergies have the 
following structure 



E'n = '-^[(^- ^sY - , (13) 

^» = ^ — — 

where w = a/1 + Q\ = 1/ cos 7, and 5 = ±1 stand for up and down spins. We emphasize 
that hereafter, the terms up and down (labeled respectively | and |) refer to directions in 
the local frame, as illustrated in Fig. 1. 
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FIG. 1: (Color online) (a) Schematic illustration of the spin orientation exhibited by the eigenstates 
of an ideal one-dimension ring, (b) and (c) show schematically the energy spectrums for the total 
angular momentum quantum number k at different spin-orbit coupling strengths w = cos~^7. 

1. The weak SOI limit 

In the weak SOI limit, i.e. if (7 — > 0) and in absence of the static external flux the eigenen- 
ergies attain the forms Ej^ = ^n"^ and -E^ = + 1)^ which seems inconsistent with 

the spin-degenerate ground-state formula En = (associated with Hq = 2m*o? (^"^y)^)- 

The resolution of this apparent inconsistency is as follows, -ffsoi commutes in a nontrivial 
way with K = ^ Sz which is the z component of the total angular momentum, and 
with S^^ = Sx sin 7 cos ip + Sy sin 7 sin + Sz cos 7 which is the spin component in the di- 
rection determined by the angles 7 and (p. Also we can show that [K, 5*^^] = 0. The 
simultaneous eigenfunctions of H, K and S"^^ are the function given by Eq. ([8]). To ro- 
tate the quantization axis of 5*^ to the direction 5*^^, a SU(2) transformation is necessary, 
i.e. US-JJ~^ = S^^, where U is the transformation matrix, Uu = cos'y/2exp{—iip/2), 
U12 = — sin7/2 exp(— 'iy9/2), U21 = sin7/2 exp(iy9/2), and U22 = cosj/2exp{i(p/2). Under 
the transformation of the wave functions of SOI become \E'' = f/~^\E'f (</?) = e*^""^^/^'''^X5) 
where xs = (Ij 0)"^ stands for spin-up state and (0, 1)^ for spin-down state. The exponential 
factor n+1/2 is the quantum number for the z component of the total angular momentum, 
referred to as K\l/^(<y9) = k\1/^(<^), where k = n + 1/2. n is fixed in the process of the SU(2) 
rotation. On the other hand, the wave function for Hq, the Hamiltonian for a ring without 
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SOI, is usually set as \l/o = e*"'^X5-~ Comparing and \l/o, it is clear that the exponential 
factors in and \[^o are different, and n in \l/o is the quantum number of orbital momen- 
tum. This difference manifests itself in the limit 7 = (i.e. SOI is zero, -ffsoi ~^ Hq): the 
wavefunctions "^^iv) ^'^ S*^ '^^^^ i'^ ^0 because they are the eigenfunctions for different 
sets of commuting observables. The simultaneous set for \l/o is H, and Sz- 

Despite this explanation we may still wonder what is the physical effect of splitting of 
the energies on n axis^^ in the limit of vanishing SOI. We inspect therefore the persistent 
charge current (CC) caused by SOI and take the limit of absence of SOI in the last step. 
The partial charge persistent current due to the state labeled by n and S in the presence of 
SOI and a static external magnetic flux is 

T f 1 — Sw\ ,^ 
InS = -e<^/o ( ^ + 0^ + — 2 — ) ' ^ ^ 

where Iq = 2£'oa/0o is the unit of current. When = 0, the CC due to the particle in 
the n level is Ins = ~n for spin up and Ins = —{n + 1) in the limit of zero SOI (we 
ignored e,^ and scaled the current by Iq). The total CC is Is = Ylm-^ns- For a distribution 
of up spins, the occupied states are n = 0, ±1,±2,---; for spin-down particles, they are 
n = — 1, (0, — 2), (1, — 3), ■ ■ ■ ; the numbers in the same parenthesis "(...)" indicate states 
with the same energy. Hence, as expected, we conclude that the total CCs for up and down 
spins carriers are zero in the limit of vanishing external magnetic field and SOI. Note, that 
the lowest occupied state for down spin is not at n = 0, but at n = —1. This means the 
velocity of the particle with down spin at = — 1 is zero; in fact n + 1 stands for the 
quantized velocity for spin-down particles. Thus, shifting the energy spectrum in the down 
spin channel to the right by 1 quantum number, the energy spectrums for different spins 
become the same on the velocity axis which is consistent with the physical picture^. 
In presence of SOI the energy spectrum is rewritable as 

E'^ = ^[(K - Sw/2 - 0/0o)^ - gy4] = ^[(. - i - - Qm (16) 

For = the energies associated with a total angular momentum quantum number k are 
spin-split; up spins have lower energy than down spins as illustrated in Fig. 1(b) and 1(c). 
For the case of zero SOI, i.e. for w = 1 (cf. Fig. 1(b)) electrons with spin-up and k, are 
degenerate with electrons having spin-down and k — 1 quantum numbers; the energy is given 
hj K — 1/2, meaning that the energy levels are spin-degenerate. 



B. Time-dependent wave functions and energies after pulse irradiation 



As have been shown in detail^'i^'^'^, upon applying at t = a half-cycle pulse with a 
duration r^, the time-dependent electronic states of the ring develop as 

(y,, to = T,) = vl/f (y,, t = 0)e-™^'^. (17) 

The parameter 

a = eap/h,p=- E{t)dt (18) 







is the action (in units of action) taken over by the carriers from the pulse EM field. The 
range of validity of the the solution (ITTl) has been discussed in Refs. jlGlllTllssI (for a general 



discussion of the properties of the time development operator we refer to the work |37l|). 
The coherent state (fT7|) is not an eigenstate of the ring for t > to, in fact a state initially 
labeled by the quantum numbers no and 5*0 is expressible in terms of the ring stationary 
eigenstates as 

K>,t) = -^$^Cf(no,^o,t)e^("-^V2Ve-^i^„^V^|z.^), (19) 



where Ip^) = h'^{'j,ip). From Eq. (fTTI) we deduce the expansion coefficients 



^SSo^nno, for t < 0, 



(20) 



SsSoi""" "Jno-n(")' fo^^ t > to, 

The energy at time t is then given by the relations 

= (v&S(^,t)|/^|vi/5(¥.,t)) = ^M^S(¥',i)i^i^S(¥',^)) 

= 5]£;„"|Cf(no,^o,t)p. (21) 

ns 

Substituting Eq. ([13]) in Eq. we find 

E^m = ^ {(„„ - .s.r + ; (22) 



> („) = < 0) + (23) 

Here 6(t) is the Heaviside step function and 0:5 is given by Eq. f|T4l) . 
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IV. DIPOLE MOMENT GENERATED BY THE PULSE UNDER EXTERNAL 
STATIC MAGNETIC FIELD 



Having identified the time-dependent spectrum and eigenf unctions, we focus now on tlie 
cliarge-dynamics and tlie induced polarization. To tliis end we inspect tlie cliarge localization 
parameter, defined slsf^^^^^ 

(cosy.)g(t) = / dv;|^S(yp,t)|2cos¥., 
Jo 

ns 

From Eq. (1131) and relation (I20p for the coefficients, the localization parameter is deduced 
after some algebra to be 

(cos(^)f»(t) =ah{n) sin(6) cos[2(no - xsM (25) 
h{Q) = Jq{Q) + J^iQ), b = uJot/2, and Q = a^2{l - cos(26)). 

Here J„ is a Bessel function with index n. The partial photo-induced dipole moment /i^°(t) 
associated with the initial state with quantum numbers no, So and the total HCP-induced 
dipole moment along the x axis fi^° (t) for the initial spin 5*0 read 

/xS(t) = ea(cos^)f°(t), (26) 
^^''it) = E /K,5o,iV,t)/xg(t). (27) 

f stands for the nonequilibrium distribution function whose derivation requires the solution 
of the kinetic equations. In principle we can employ our recent approach based on the density 
matrix'20 but more detailed knowledge on the spin-dependent decay channels in confined 
geometry is needed. Here we inspect the zero temperature behaviour of the induced dipole 
moment and the associated emission spectrum. We expect the quahtative features of these 
physical quantities to persist at finite temperatures, as we demonstrated for the case of 
vanishing SOI—. 

A. Spectral analysis 

In the following we focus on the spectral properties. The spin-orbit interaction breaks 
the energy-degeneracy of n and — n states. As evident from Eqs. (I14ll22p the spectrum, 
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posses a symmetry axis (SA) located at xs = 0/0o — (1 — Sw)/2 (the global additional 
pulse-associated energy and SOI-induced energy shift do not affect this symmetry), i.e. the 
static magnetic flux and the SOI act as an effective magnetic field. In Eq. fl22p uq is an 
integer, hence it is advantageous to introduce the integer parameter Is as the nearest integer 
that is less than xs and define (cr =1 or |) 



A^ = xs- h 



(28) 



whose meaning is illustrated in Fig. 2(c) and 2(d). We distinguish four cases (a), (b), (c) 
and (d) corresponding to A^. = 0,1/2, A^. G [0,1/2] and A^. G [1/2,1], respectively (cf. 
Fig. 2). Act = 1 is equivalent to Act = 0, thus Act is periodic with changing SA and varies 
within the fundamental interval [0, 1]. Furthermore, we introduce 



A. 



|l/2- ActI 



(29) 



as the distance between the SA and the half integer axis. These four cases in Fig. 2 are 
valid for up and spin-down states. 



(50 

i-H 

u 




FIG. 2: (Color online) (a) SA has an integer value; (b) SA is a half integer; (c) SA lies in the region 
(/, / + 1/2); (d) SA is in the region (l + 1/2, / + 1). The red dash lines are the SA in each case. 
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1. Spinless particles 



The two spin states, up and down, need to be considered. Allowing each of these two 
spins to occupy the four configurations shown in Fig. 2, results in 16 combinations. For 
simplicity, we consider at first only one kind of spins and states which can be any of the 
four cases depicted in Fig. 2. Depending on the total number of electrons in the ring N two 
situations are distinguished: 

(1) N is an even integer. For case (a) in Fig. 2 (for brevity we refer hereafter to 
cases (a), (b), ■ ■■; and write Is = I, A^- = A, and xs = x) the occupied states are at 
Uq = |(/ — m), / — (m — 1), ■ ■ ■ , / — 1, /, / + 1, / + 2, ■ ■ ■ , / + (m — 1), |(/ + m). Here, e.g. 
|(/ — m) means half occupation on the state characterized by the quantum number / — m. 
The dipole moment in this case /i^^^ ( A^, t) reads 

lj,l^){N,t) = ae{t - to) h{n) sin cos[2{no - x)b], 

no 

= ae{t-to)h{n) sin{Nb) cos(6). (30) 

For (b), (c) and (d) cases, we have uq = / — (m — 1), ■ ■ ■ ,1 — 1,1, l+l, 1+2, ■ ■ ■ , /+(m — 1), l+m. 
Accordingly we deduce 

fxlf,^^){N,t) = ae{t - to)h{n) sin{Nb) cos[(l - 2A)6], (31) 

where A = x — 1. In (a) case, A = 0, then Eq. fl3T]) will reduce to Eq. fl30|) . The dipole 
moment for even occupation has in general the form 

^i'{N, A, t) = ae{t - to)h{n) sin{Nb) cos[(l - 2A)6]. (32) 

(2) A^ is an odd integer. Similar steps as in the proceeding case can be performed, leading 
us to conclude that 

/i°(A^,A,t) = ae{t - to)h{n) sin{Nb) cos[(l - 2A)6]. (33) 
2. particles and (/) 

For spin-i particles the energy spectrum is spin-dependent. The position of the spectrum 
SA is different for the two spin states, i.e. X| = (p / (po + {w — 1) /2 and xi = x^ — w. Thus, the 
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relative distance between these two SAs depends on the external flux and SOI parameter 
w. Scanning these parameters the spectrums of up and down spins are tuned to any of the 
four cases shown in Fig. 2. In what follows we use for cases in Fig. 2 that correspond to 
the spin-down and the spin- up states the symbols (ij) (i,j=a,b,c,d); e.g. the case (ab) means 
the down-spin spectrum corresponds to the case (a) whereas the up-spin spectrum is as in 
case (b); hence there are 16 combinations of such pairs to be considered. In the following, 
we consider four different occupations in the ring for up and spin-down states and these 16 
combinations in detail. 

(1) Even number of pairs. For N — Am, where m is an integer {N is the total number of 
electrons in the ring) there are 2m spin-up particles and 2m spin-down particles occupying 
the respective spectrums, meaning that 

l,''{N,t)=i,%N/2,A,,t). (34) 

(2) Odd number of pairs. UN — Am + 2 then 2m -|- 1 spin-up particles and 2m -|- 1 
spin-down particles populate the respective spectrums and 

/x^(7V,t)=/x°(7V/2,A,,t). (35) 

(3) Even number of pairs plus an extra particle. Here we write N — Am + 1. i) For 2m 
spin-up particles and 2m -|- 1 spin-down particles we find ii^{2m,t) — //^(2m, Af, f), and 
//■'•(2m + l,t) — //°(2m -|- 1, Aj^,i). Cases (ab), (ac), (ad), (cb), (db) belong to this type, 
ii) If there are 2m -|- 1 spin-up particles and 2m spin-down particles in the ring we infer 
fj,^ {2m + l,t) = fi°{2m + l,A|,t), and /x^(2m, t) = //'^(2m, Aj^, t). Cases (ba), (ca), (da), 
(be), (bd) belong to this category, iii) If we have 2m spin-up particles and 2m spin-down 
particles plus one extra particle we shall analyze the populated state of the extra particle. 
E.g., cases (aa), (bb), (cc), (dd), (cd) and (dc)are possible situations. Careful calculation of 
all cases results in the formula 

pi^^ = ae{t - to)h{n) sin(6) cos[(2m + 1 - 2A^)6], (36) 

where Ao- > A^. 

(4) Odd number of pairs plus an extra particle. For N — Am -|- 3 we distinguish: i) 
For 2m -|- 1 spin-up particles and 2m -|- 2 spin-down particles we infer ^^{2m + l,i) = 
//°(2m -M, A|,i), and /x^(2m-F2,t) = //^(2m 2, A^, i). Cases (ab), (ac), (ad), (cb) and 
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(db) belong to this type, ii) For 2m + 2 spin-up and 2m + 1 spin-down particles we deduce 
/iT(2m + 2,t) = fi%2m + 2, A|,t), and /i^(2m + l,t) = fi°{2m + 1, A|,t). Cases (ba), (be), 
(bd), (ca) and (da) are examples for this situation, iii) For 2m -|- 1 spin- up and for 2m -|- 1 
spin-down particles plus an extra particle (cf. (bb), (cc), (dd), (cd) and (dc)) we find for 
the dipole moment 

/ig^ = aQ{t - to)h{n) sin(6) cos[(2m + 1 + 2A^)b] (37) 

where < Ag-. 

3. spin ^ particles with ((p = 0) 

If = then = X| = — and A| -|- A| = 1 applies. The symmetry axes for up 
and down-spin states have the same distance from the nearest integer axes to the left and 
to the right sides to the SAs. The dipole moments are spin degenerate. 

(1) Even number of pairs and = 0. For = 4m we find fi\N, t)=/i|(2m, A|, t). 

(2) Odd number of pairs and = 0. For = 4m-|-2 we conclude fi\N, t)=fi°{2m + l, A|, t). 

(3) Even number of pairs plus an extra particle. If N = 4m -|- 1, only possible combinations 

), (bb), (cd) and (dc). The spin degeneracy of the extra particle is caused by the 
crossing of the energy levels with opposite spins. The dipole moment reads 

^T(7v, t) = ^aQ{t - to)h{n) {sin(2m6) cos[(l - 2A^)b] + sin[(2m + 1)6] cos(2A|6)} . (38) 

If A I = 0, Eq. (l38l) delivers the dipole moment for the ) whereas A-f = 1/2 applies 

for the case (bb) and others for case (cd) and (dc). 

(4) Odd number of pairs plus an extra particle. For A^ = 4m -|- 3, only (aa), (bb), (cd) and 

(dc) are applicable. The dipole moment is 

^\N, t) = ia0(t - to)h{Q) {sin[(2m + 1)6] cos[(l - 2A|)6] + sin[(2m + 2)6] cos(2A|6)} . 

(39) 

In the event A| = 0, Eq. (l39l) yields the dipole moment for case (bb) whereas A| = 1/2 is 
valid for the case (aa) and others for case (cd) and (dc). 
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V. NUMERICAL RESULTS AND DISCUSSIONS 



A. Experimental feasibility and general remarks 

In this section we present and analyze numerical results for the HCP-induced polarization 
of a ballistic quantum ring. In view of an experimental realization it is important to identify 
the realistic range of the parameters such as the strength of the SOI and the associated 
quantities. The range of the ring size and external field strength are chosen according to 
current experimental feasibility. 

The Rashba SOI was already investigated for numerous semiconductor quantum wells 
such as In2;Gai_a;As/InP quantum wells^°, Ino.53Gao.47As/Ino.52Gao.48As heterostructures^, 
and GaSb/InAs/GaSb quantum wells^. Spin-interference effects in a ring with Rashba 



SOI which was built in a InGaAs/InAlAs heterostructure was the subject of Ref. 43| . 
and the Aharonov-Casher phase was inspected in II- VI semiconductor quantum rings, such 
as HgTe/HgCdTe ring^. In the context of the present work it is important to estimate 
the realistic range for the spin orbit angle 7 for the relevant semiconductor materials: For 
Inx.Gai_x'As/InP quantum well^, varies in the range [0.7 x 10~^^el^-m, 1.1 x lO^^^eV -m] 
when an applied gate voltage varies from +1.5 V to -2.5 V; this corresponds to 7 be- 
ing in the range [-34°, -47°], [-74°, -79°], and [-82°, -85°] for rings with 100 nm, 500 
nm, and 1 /im radius respectively (m* = 0.037mo, mo is the free carrier mass). For 
Ino.53Gao.47As/Ino.52Gao.48As materials"^^ , varies in [0.5 x 10~"'^^el^ ■ m, 1.0 x lO^^-'^eV^ ■ m] 
if the external gate voltage changes from +1.5 V to -1 V; correspondingly 7 varies in 
[-33°, -52°], [-73°, -81°], and [-81°, -86°] for rings with 100 nm, 500 nm, and 1 /xm 
radius respectively (m* = 0.05mo). For GaSb/InAs/GaSb quantum well^ the values of 7 
are on the same order as above. 

With regard to available pulses, a wide range of pulse durations and strengths have been 
realized^'2^'^'2^'^. As clear from the above analysis the pulse parameter which is decisive 
for the electron dynamics is a as given by Eq. ( JTSll . For pulse with a sin-square shape 
and a duration of r^^ = 1 ps we achieve for a ring with radius 1 yum a transferred action of 
ai = 0.1, a2 = 1, and 0:3 = 10 if the peak electric field strength is tuned to respectively 
El = 1.32 V/cm, E2 = 13.2 V/cm, and = 132 V/cm. It is this range of a which we use 
in the present numerical calculations. In the figures below we just provide a. 
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From a general point of view we can expect four frequency scales to be relevant for the 
time-dependent charge and spin dynamics: 

1) . The global energy scale is set by the size of the system. Hence, the fundamental frequency 
is given by c^o = h/{m*a^) and the natural time scale is tp = Att/ujq which for the ring sizes 
at hand is tens of picoseconds. 

2) . As for any fermionic system, the Fermi energy Ep sets the scale for the fast (femtosecond) 
charge dynamics associated with excitation near Ep- Since we are dealing with an isolated 
ring the Fermi energy is expressible in terms of the number of particles and the relevant 
frequency is therefore Nujq. 

3) . A Further frequency ur is associated with SOI-influenced dynamics; uji = 2afi/{ha). 
For our systems ujr can be on the order of ujq opening the possibility for controlling quantum 
interferences, e.g. by tuning the strength of SOI (via a gate voltage). 

4) . Further modifications are brought about by the applied pulse which (for the pulses of 
interest here) induces a multitude of excitations near Ep with further associated frequencies, 
as detailed below. 

These physical expectations are confirmed by the general structure of the calculations for 
the induced dipole moments according to Eqs. fl25|30|32ll33|36|37l) 




B. Zero static magnetic field 

At first we investigate the case of vanishing static flux (0 = 0). In Fig. 3 the time- 
dependence of the dipole moment is shown for different spin-orbit angles 7 and pulse 
strengths a. Because of the spin-degeneracy we consider only one spin channel. The total 
number of particles is A^ = 4m. The time is measured in units of the system's time scale 

tp = 4:71171*0^/11. 

The four time scales mentioned above are clearly visible. As inferred from Eqs. ( l30l 
[32l [33l [36l [37j) the net time- dependent dipole moment grows linearly with the pulse field 
strength at small a (note for a^l/2, we have h{Q) ^ 1). The fast oscillation in 

the dipole moment are related the transition between the levels near Ep (Rabi flopping). 
With increasing a more levels are excited giving rise to higher contributing harmonics (cf. 
a = 1 and a = 3 in Fig. 3 and also Fig. 4). 

The SOI strength (quantified by 7) has a dramatic influence on the low frequency mod- 
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FIG. 3: (color online) The time-dependence (in units of tp = Att/ujq) of the net pulse-induced dipole 
moment fi{t) (in units of ea) in the up or down-spin channel at different 7 and a as depicted on 
the figure. The total number of particles is = 100 and ^ = 0. 

ulation of the dipole-moment envelope. The low frequency is associated with the difference 
between the frequencies of the involved levels near Ep, which is on the scale of ujq. In fact, as 
inferred from Eq. (132|) 7 can be tuned as to influence the phases of the involved wave func- 
tions changing thus the interference pattern and removing eventually the slow oscillations 
altogether. This happens at 7 = —60° as shown in Fig. 3. When I7I is increased further 
the slow modes appear again. This is insofar important as 7 can be modified by an external 
gate voltage offering thus the possibility of engineering the emission spectrum via an applied 
static electric field and opening the way for testing experimentally our theoretical predic- 
tions. For this reason we inspect the power spectrum^^ produced by the non-equilibrium 
charge oscillations in the QR by evaluating 

/oo 
li{t)e-'^Ut. (40) 
-00 

Fig. 4 shows the power spectrum for different strengths of the pulse a. As evident from 
these calculations the frequency scale is set by tUo/S. In Fig. 4 the number of particles 
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Frequency (co/((»y2)) Frequency ((£>l{a>J2y) 

FIG. 4: (color online) The emission spectrum for different pulse strengths (quantified by a) and 
SOI strengths (indicated by 7). In (a) a is varied at 7 = whereas in (b) 7 is variable at fixed 
a = 0.01; in both cases N = 100. 

occupying the single spin states is N=50, for small a only the state at Ep is excited leading 
to the appearance of two frequencies at (A^ — l)uJo/2 and (A^ + 1)ujq/2. With increasing a, 
more levels are excited and correspondingly further frequencies in unit of ujq/2 emerge at 
47, 45, ■ ■ ■ and 53, 55, 57, ■ ■ ■ (see Fig. 4(a)). In this context we note that a very short HCP 
contains almost all frequencies. Nevertheless, at low intensities only a limited number of 
states in the ring can be excited. The reason is obvious from Eq. ( l23l) . The highest energy 
level achieved upon excitation is Ep + hw^a^ Hence there is an excited energy cut-off set 
by the field intensity (available photons) and limits consequently for a certain a the number 
of possible frequencies observable in the power spectrum (as seen in Fig. 4). 

Fig. 4(b) shows how the SOI shifts of the frequencies: with increasing I7I the frequency at 
49 moves towards a higher frequency while the frequency at 51 moves to a lower frequency. 
For the angle 7 = —60° we infer A = 1/2 and the two frequencies merge into one frequency. 
Further increasing the SOI strength, i.e. A > 1/2, the frequency peak from 51 continues 
moving to 49, while the frequency peak from 49 approaches 51 until they coincide for A = 1 
(or A = 0). This behavior of frequencies is repeated periodically with increasing SOI. 
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FIG. 5: (color online) The time-dependence of dipole moment is shown at different occupation 
numbers. Here 7 = —70°, a = 0.1 and m = 25. Units are as in Fig. 3. 

As detailed above the dipole moment depends sensitively on the occupation numbers. 
Fig. 5 shows an example of the dipole-moment dynamics for = Am + 2 contrasted with 
the case N = Am. We remark that since A + A = 1/2 we have A = 0(i) for A = i(0). Thus 
for = Am, when the dipole moment increases with A in A G [0, |] the dipole moment 
decreases for the case A^ = Am + 2. Furthermore, the parity of the occupation number is 
very important for the property of dipole moment of the ring. E.g., for A ^ A, the case 
A^ = Am {N = Am + 1) behaves similarly to A^ = Am + 2 (A^ = Am + 3) in which cases the 
occupations are even (odd). However, the even case is qualitatively different from the odd 
case (not shown in Fig. 5). If A 7^ A the oscillations of the dipole moment are different in 
all the cases depicted in Fig. 5. 

C. Finite static magnetic field 

A static external magnetic field lifts the spin degeneracy of the dipole moment. The 
induced dipole moment dynamics becomes spin dependent (cf. Fig. 6). The slow- frequency 
shift between the dipole-moment oscillations in the up-spin and down-spin channels is readily 
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FIG. 6: (color online) The time-dependence of the dipole moment for up and down spins under a 
nonzero static magnetic field are shown in (a) and (b). The power spectra corresponding to the 
up and down spins are shown in (c) and (d) respectively. Here (j) = 0.25, 7 = —70°, a = 0.1 and 
N = 100. 

understood from the energy splitting (cf. Eq. ( l23l) ). The power spectrums for the up and 
down-spins dipole oscillations (Fig. 6(c) and 6(d)) reveals clearly the frequency shift. 

For more insight into the role of the SOI we investigate the spin-resolved local charge 
density before and after the pulse. The probability density associated with the level labeled 
by no and 5*0 is p^°{ip,t) = \'^^°{{p,t)\'^. Before the pulse the charge density is a unit charge 
uniformly distributed around the ring, i.e. po = Upon applying the pulse we find 

pg(</.,t>0)=po+ Apg(^,t>0), (41) 

where the second term is the spin-resolved, field-induced charge density variation (ICDV) 

which we evaluated and find 

^ 00 

Apf"(^,t>0) = -{y^J2n+i[2asm{2n + l)b]cos{2n+l)[if-2{no-Xso)b] 

n=0 
00 

+ ^J2n['2asm{2nb)] cos2n[ip - 2{no - X5J6]}. (42) 

n=l 

The total spin-resolved charge density is obtained from a sum over the occupied levels 
p^°{(p,t > 0) = NsoPo+ A p^°{ip,t > 0), where A^s^ is the total particle number in the 5*0 
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spin channel. The term A p^''{ip,t > 0) for the even pair occupation case is 

oo ^ 

A p^°{<f,t > 0) = -{V 7— — [J2„(^^in) + ^2n+2(^^in)] sinlNs.i^n + 1)6] cos(2n + l)^Soit) 
vr ^ — ^ Zn + i 

n=0 
oo ^ 

+ J^^[-^2n-i(^2n) + J2n+i{^2n)] sm{Nso2nb) COS 2rup Soit)} , (43) 

n=l 

where = aasin(2?2 + 1)6, f22„ = aa;sin(2n6), and '^Soit) = ~ (1 ~ 2A5q)6. 

90 (l/2)Ap (y=0°, <t,=0) 
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FIG. 7: (color online) A snap-shot of the pulse-induced charge density variation (ICDV) as a 
function of the azimuthal angle if at the time t/tp = 2.005. The pulse parameter is a = 0.1, and 
N = 100. Full (open) dots show the ICDV in the up (down) spin channel. Solid curve stands for 
the spin-averaged ICDV. 

Eq. ( H3!) indicates a spin-dependent phase of A caused by the interplay of SOI, the 
static magnetic field, and the pulse field. Interestingly the phase shift evolves with time. 
Fig. 7 shows the time-evolution of the spin-dependent ICDV. In all cases depicted in Fig. 
7 the static flux is absent. For 7 = we observe, as expected, how the pulse kicks the 
charge density along the pulse-polarization axis. For the given time moment, the missing 
charge density around (f = it is pushed to the region around ip = 0. For a finite spin orbit 
interaction (I7I > 0) we observe a spin splitting of ICDV, meaning that the pulse induces 
temporally and locally a finite spin polarization P = p-^ — pi, even though the system 
is initially paramagnetic. The time integral of the pulse-induced, local spin polarization 
vanishes however. Technically, we infer that the SOI results in a SU(2) flux that produces 
opposite phase shifts on the azimuthal angle with the same magnitude for a zero-static 
magnetic field, i.e. ^p^^) = =f (1 — 2A^)b. These shifts cause a rotation around the ring 
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of the symmetry axes of the spin-up and the spin-down densities respectively clockwise and 
anticlockwise (see Fig. 7). When 7 = —60° (A| = 1/2), the up and down ICDV merge into 
one curve after one period. The periodic rotation is subject to the condition A| = 1/2. The 
spin-averaged ICDV, i.e. A p = A p|+ A p|, is always symmetric with respect to x axis. At 
(y9 = and 7r ICDV is spin-degenerate for all times (-P(t) = 0). 

The time evolution of the spin-resolved ICDV is depicted in Fig. 8. The symmetry 
axes of the pulse-induced spin-up and spin-down ICDV rotate around the ring respectively 
anticlockwise and clockwise with time (which is the opposite behaviour when increasing 
SOI). Accordingly the total ICDV oscillates along the x axis and possesses a left-right 
symmetry with respect x. The local and temporal spin polarization P is symmetric to y 
axis and oscillates along it with the same frequency as p| (or p|). 




FIG. 8: (color online) The time and spatial dependence of the ICDV shown in polar coordinates 
in (a), (b), (c) and (d) at different times after the pulse. Filled circles represent the up spin ICDV 
and the open circles indicate the down spin ICDV. The parameters are a = 0.1, = 0, 7 = —40° 
and N = 100. 
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VI. CONCLUSIOS 



In summary, we investigated the dynamical response of a quantum ring with spin-orbit 
interaction upon the apphcation of a hnearly polarized time-asymmetric weak electromag- 
netic pulse. It is found that the dipole moment along the pulse polarization axis is spin 
degenerate when there is no external magnetic field. The dipole moment oscillates with 
time and the SOI provides an envelope function or shifts of the oscillation frequencies of 
dipole moment. Stronger pulse fields can excite higher and lower harmonics. And the enve- 
lope functions for different parity of the occupation number on the ring are different. When 
a static external magnetic field is applied, the spin degeneracy is removed. The spatial and 
temporal dependence of the pulse induced charge-density variation (ICDV) indicates that 
the SOI results in a SU(2) fiux leading to a sphtting of the phases of the up and down spins 
states. The symmetric axes of the ICDV for the up and down spins arc rotated equally in 
clockwise and anticlockwise when increasing SOI. The total ICDV and the local and tem- 
poral polarization of the charge density are symmetric to respectively the light polarization 
axis and the axis perpendicular. The pulse-induced polarization is experimentally accessible 
by measuring the power spectrum of the emitted radiation. 

We thank A.S. Moskalenko and A. Matos-Abiague for fruitful discussions. The work is 
support by the cluster of excellence " Nanostructured Materials" of the state Saxony- Anhalt. 
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